Induced fermionic current in toroidally compactified spacetimes 
with applications to cylindrical and toroidal nanotubes 



S. Bellucci 1 * A. A. Saharian 2 f V. M. Bardeghyan 2 

1 INFN, Laboratori Nazionali di Frascati, 
Via Enrico Fermi 40,00044 Frascati, Italy 

2 Department of Physics, Yerevan State University, 
1 Alex Manoogian Street, 0025 Yerevan, Armenia 

August 24, 2010 



Abstract 

The vacuum expectation value of fermionic current is evaluated for a massive spinor 
field in spacetimes with arbitrary number of toroidally compactified spatial dimensions in 
the presence of a constant gauge field. By using the Abel-Plana type summation formula 
and the zeta function technique we present the fermionic current in two different forms. 
Non-trivial topology of the background spacetime leads to the Aharonov-Bohm effect for 
the fermionic current induced by the gauge field. The current is a periodic function of the 
magnetic flux with the period equal to the flux quantum. In the absence of gauge field it 
vanishes for special cases of untwisted and twisted fields. Applications of general formulae 
to Kaluza-Klein type models and to cylindrical and toroidal carbon nanotubes are given. In 
the absence of magnetic flux the total fermionic current in carbon nanotubes vanishes, due 
to the cancellation of contributions from two different sub-lattices of the hexagonal lattice of 
graphene. 



PACS numbers: 03.70. +k, ll.10.Kk, 61.46.Fg 



1 Introduction 



In many physical problems we need to consider some model on a background of manifold having 
compact spatial dimensions along which dynamical variables satisfy some prescribed periodic- 
ity conditions. Incomplete list of applications where the topological effects play an important 
role includes Kaluza-Klein type models, supergravity and superstring theories. From an infla- 
tionary point of view, universes with compact dimensions, under certain conditions, should be 
considered as a general rule rather than an exception [lj. Models of a compact universe with 
non-trivial topology may play an important role by providing proper initial conditions for in- 
flation. An interesting application of the field theoretical models with non-trivial topology of 
spatial dimensions appeared in nanophysics recently [2] . The long- wavelength description of the 
electronic states in graphene can be formulated in terms of Dirac-like theory in 3-dimensional 
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spacetime with the Fermi velocity playing the role of a speed of light (see, e.g., Refs. OS]). 
Single- walled carbon nanotubes are generated by rolling up a graphene sheet to form a cylinder 
and the background spacetime for the corresponding Dirac-like theory has a topology R 2 x S 1 . 
The compactification in the direction along the cylinder axis gives another class of graphene 
structures called toroidal carbon nanotubes with the background topology R 1 x (S 1 ) 2 [5]. 

The compactification of spatial dimensions leads to a number of interesting field theoretical 
effects which include instabilities in interacting field theories, topological mass generation and 
symmetry breaking. In quantum field theory the boundary conditions imposed on fields along 
compact dimensions change the spectrum of vacuum fluctuations. The resulting energies and 
stresses are known as the topological Casimir effect. (For the topological Casimir effect and its 
role in cosmology see [6]- [10] and references therein.) Note that Casimir forces between material 
boundaries are presently attracting much experimental attention [TT]. In the Kaluza-Klein type 
models this effect has been used as a stabilization mechanism for moduli fields which parametrize 
the size and the shape of the extra dimensions. The Casimir energy can also serve as a model 
of dark energy needed for the explanation of the present accelerated expansion of the universe 
(see |12j and references therein). 

The effects of the toroidal compactification of spatial dimensions on the properties of quan- 
tum vacuum for various spin fields have been discussed by several authors (see, for instance, 
[6]-|13] and references therein). One- loop quantum effects for the scalar and fermionic fields in 
de Sitter spacetime with toroidally compactified dimensions are studied in Refs. [HI [15]. I n P re " 
vious papers [Ml [T7] we have investigated the fermionic condensate and the vacuum expectation 
value of the energy-momentum tensor for a massive spinor field in higher-dimensional space- 
times with toroidally compactified spatial dimensions. These expectation values are among the 
most important quantities that characterize the properties of the quantum vacuum. Another 
important characteristic is the vacuum expectation value of fermionic current. Although the 
corresponding operator is local, due to the global nature of the vacuum, this quantity carries 
an important information about the global properties of the background spacetime. In addi- 
tion to describing physical structure of the quantum field at a given point, the current acts 
as the source in the Maxwell equations. It therefore plays an important role in modelling a 
self-consistent dynamics involving the electromagnetic field. 

In the present paper, we investigate one-loop quantum effects on the fermionic current arising 
from vacuum fluctuations of a massive fermionic field on the background of spacetimes with an 
arbitrary number of toroidally compactified spatial dimensions. We will assume generalized 
periodicity conditions along the compactified dimensions with arbitrary phases and the presence 
of a constant gauge field. A non-zero gauge field defined on topologically non-trivial background 
leads to the Aharonov-Bohm effect for the vacuum expectation value of fermionic current. Note 
that fermionic current in spacetime with non-trivial topology induced by a cosmic string has 
been investigated in Refs. [18] . 

This paper is organized as follows. In the next section, we consider vacuum expectation value 
of the fermionic current in the background spacetime with spatial topology R p x (S 1 ) q in the 
presence of a constant gauge field. The corresponding expression is derived by using the Abel- 
Plana type summation formula. An equivalent representation is obtained in Section [3] within 
the framework of the generalized zeta function approach. In Section 0] we apply the general 
formula for the evaluation of the fermionic current in cylindrical and toroidal nanotubes within 
the framework of Dirac-like model for electrons in graphene. Main results are summarized in 
Section [5j 
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2 Vacuum expectation value of the fermionic current 



We consider the quantum fermionic field tp on a background of (D+l)-dimensional flat spacetime 
with spatial topology R p x (S 1 ) 9 , p+q = D. The Cartesian coordinates along uncompactified and 
compactified dimensions are denoted ) and z q = (z p+1 , . . . , z D ), respectively. 

The length of the Z-th compact dimension we denote as L/. Hence, for coordinates one has 
— oo < z < oo for I = 1, . . . ,p, and ^ z l ^ L\ for I = p + 1, . . . , D. We assume that along the 
compact dimensions the field obeys the generic quasiperiodic boundary conditions, 

i/>(t, z p ,z q + L iei ) = e 2lTiai iP(t, z p , z q ), (1) 

with constant phases a/ and with being the unit vector along the direction of the coordinate 
z , I = p + 1,. . . ,D. Condition ([I]) includes the periodicity conditions for both untwisted and 
twisted fermionic fields as special cases with a/ = and a\ = 1/2, respectively. As is discussed 
below, the special cases a\ = 0, ±1/3 are realized in nanotubes. 

Dynamics of the massive spinor field is governed by the Dirac equation 

i^D^-mip = , D fl = d fl + ieA fl , (2) 

where is the vector potential for the external electromagnetic field. In the discussion below we 
assume that = const. Though the corresponding magnetic field strength vanishes, the non- 
trivial topology of the background spacetime leads to Aharonov-Bohm-like effects for physical 
observables. In particular, as it is shown below, the expectation value of fermionic current 
depends on A^. In the (D + l)-dimensional spacetime, the Dirac matrices are N x N matrices 
with N = 2^ D+1 ^ 2 \ where the square brackets mean the integer part of the enclosed expression. 
We assume that these matrices are given in the Dirac representation: 

7°=(J -0'^(-°j '•')•« " 

From the anticommutation relations for the Dirac matrices one has cr^a^ +ov<7+ = 25 In the 
case D = 2 we have N = 2 and the Dirac matrices are taken in the form 7 M = (<7p3, zcpi, i<rp2); 
with ap^ being the 2x2 Pauli matrices. We are interested in the effects of non-trivial topology 
on the vacuum expectation value (VEV) of fermionic current = ipj^ip, where ip is the Dirac 
conjugated spinor. Note that the fermionic condensate and the VEV of the energy- momentum 
tensor in the model under consideration were evaluated in Ref. |16j in the absence of gauge 
field. 

By expanding the field operator in terms of annihilation and creation operators, the VEV of 
fermionic current is presented as the sum over all modes 

(0!^|0) = ^4- ) (x) 7 ^ ) (^), (4) 

where {V , ( g + '\Vy3 ^} is the complete set of positive- and negative-frequency eigenfunctions sat- 
isfying the periodicity conditions (pQ) along compact dimensions. Here, (3 is a set of quantum 
numbers specifying the solutions. The dependence of the eigenfunctions on the spacetime coor- 
dinates can be taken in the form e lk - r_ia; * ) with the wave vector k. From the Dirac equation for 
the positive- and negative-frequency solutions we find 

^ ~ A f 6 { ^^(k-eA)W +) J' (5) 
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where j3 = (k,a), = (Aq,— A), and a = (<7i, 02, . . . , <7£>)- I n these expressions w^~\ a = 
1, . . . , N/2, are one-column matrices having N/2 rows with the elements = <5<rf, an d w o- = 
ui7a • The frequency and the wave vector are connected by the relation (oj =F eAo) = (k T 
eA) 2 + m 2 for the function ip^ . 

The coefficients are found from the orthonormalization condition J d 3 xip^ + ip^ = 

80/31, where 6pp> is understood as the Dirac delta function for continuous indices and the Kro- 
necker delta for discrete ones. From this condition one finds 



A±)2 1 



{uj =F eA - mf 
(kTeA) 2 



(7) 



where V q = L p+ \ ■ ■ ■ Ld is the volume of the compact subspace. 

We decompose the wave vector into components along the uncompactified and compactified 



dimensions: k = (k p , kg), k = y k 2 + k 2 . The eigenvalues for the components along the compact 
dimensions are determined from boundary conditions (pQ): 

h = 2-K(m + ai)/L h m = 0, ±1, ±2, . . . , I = p + 1, . . . , D. (8) 

For the components along the uncompactified dimensions one has —00 < k\ < 00, I = 1, . . . ,p. 

Substituting the eigenfunctions ([6]) into the mode sum formula (jlj) and by using the properties 
of Dirac matrices one finds 

(01/10) = 4- E 1. (0) 



2V q J (2tt)p 



(0\i l \0) = — ! ^ Y , (k + eA) ' (10) 
WJ |U; 2V q J (2tt)p ^ V(k + eA) 2 +m 2 ' ^ 

with I = 1, 2, . . . , D and n q = (n p +x, . . . ,no)- In order to give a meaning to divergent expres- 
sions, it is necessary to regularize them. Here we use a Pauli-Villars gauge-invariant regular- 
ization. An alternative way is to introduce a cutoff function. Introducing regulator fields with 
large masses M s , s = 1, 2, . . . , S, for the regularized expressions one finds 

^ = ^£^ (11 > 

(Q\ji\ Q ) - A I top. V V C a Qi + eA)i_ 

(0|J " 2V q J (2n) P 2^ ^ V(k + eA)2 + M 2 ' (12) 

with Co = 1 and M = m. Under the conditions £)f=o C s M 2n = 0, n = 0, 1, . . . , [L>/2], these 
expressions are finite. After the renormalization subtractions the regulator is removed taking 
the limit M s — > 00, s ^ 1. From formula (llip it follows that the VEV of the temporal component 
of the fermionic current is renormalized to zero. Shifting the integration variable in (|12p . we 
directly see that (0|/ |0) reg = for the components with I = l,...,p. Hence, the renormalized 
VEV of the fermionic current is different from zero only for the components along the compact 
dimensions. 

Shifting the integration variables, k\ + eAi — > ki, I = 1, . . . ,p, for these components one finds 
2V q J (2n)P n?£Zq s=Q / k 2 ~ E f =p+1 [2vr(n, + a,)/^] 2 + M| 
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with r = p + 1, . . . ,D, and we have introduced the notation 



ai = ai + eAiLi/(2n). (14) 

Hence, the VEV of the fermionic current depends on components of the vector potential along 
the compact dimensions alone. We present the second term on the right of f|14[) as eAiL[/(2ir) = 
Ni + 7/, where Ni is an integer number and 7^ is the fractional part. As it is seen from formula 
(113p . only the fractional part leads to nontrivial effects. Another point to be mentioned is 
that the presence of a gauge field leads to the shift of the phases in the quasiperiodic boundary 
conditions along compact dimensions. This feature is applicable to the fermionic condensate and 
the VEV of the energy-momentum tensor as well. In particular, the formulae for these quantities 
in the presence of a constant gauge field are obtained from the corresponding formulae in Ref. 
[16] by the replacement a\ — > a\. 

The property that the VEVs depend on the phases a\ and on the vector potential components 
along compact dimensions in the combination (114ft can also be seen by the gauge transformation 

= A'^ + d fl A(x), ip(x) = if) 1 ' (x)e~ te ^ x \ with the function A(x) = A^x^. The new function 
ijj'{x) satisfies Dirac equation with A'^ = and the quasiperiodicity conditions similar to (pQ) 
with the replacement ai — > a.\. Corresponding eigenspinors are given by expressions ([5]), ([6]) 
with A = and the eigenvalues for the wave vector components along compact dimensions are 
defined by k\ = 2ir(ni + aj)/Lj. In the new gauge, the regularized VEVs are given by Eqs. (fTTj) 
and ([12]) with A = 0. The latter coincide with ([TTD and {12]) after the shift k x + eA t k u 
I = 1, . . . ,p, of the integration variables for the components along uncompactified dimensions. 

We will evaluate the VEV of fermionic current by two equivalent methods: by applying the 
Abel-Plana type summation formula and using the zeta-function technique. In the first approach 
we apply to the series over n r in Eq. (|13p the following summation formula: 



00 />0O 

V] g(n r + a r )f{\n r + a r \) = / du [g(u) + g(-u)] f(u) 



n r =— 00 



+*/ du [f(iu) - f(-iu)} ^ ^JStT - (15) 



•- l — < p *„(u+iXa r ) _ 1 ' 
" A=±l 

This formula is obtained by combining the summation formulae given in Ref. |19j (see also |16|). 
In the special case of g(x) = 1, a r = formula (|15|) reduces to the standard Abel-Plana formula 
(for the applications of the Abel-Plana formula and its generalizations in quantum field theory 
see [SHHl EE]). Taking in Eq. (fTo) 

g(x) = 27rx/L r , f(x) = \k 2 p + u 2 n r + (2irx/L r ) 2 } ^ , (16) 



with 



'q-l 



D 

E ^{m+a^/Li] 2 + m 2 , (17) 

l=p+l : l^r 

and n r q _ 1 = (n p+ i, . . . , n r _i, n r+ i, . . . ,no), we see that the first integral on the right-hand side 
of this formula vanishes. The contribution of the second integral to the regularized VEV is finite 
in the limit M s — > 00, s ^ 1, and in this term the regulator can be safely removed. By using 
the expansion l/(e y — 1) = Ylri=i e ~ ny m the integrand of the second integral, the integrals with 
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the separate terms in this expansion are evaluated explicitly and one finds 



2n(n r + 5i r ) / L 



2L r cx 



^■y^sin(27rno; r )^/kg + a;^ ~i^i(nL r y^+a^ ~), (18) 



71=1 

where K v {x) is the modified Bessel function. The equality in Eq. (jlgj) is understood in the sense 
of the renormalized value. 

By taking into account the result (j!8p . from Eq. fjl3|> . after the integration over k p , for the 
renormalized VEV one finds 

/n\-r\n\ 2NLr sin(27ma r ) p / 2 +i „ , T % , , 

v ' q n=l K r > n^_ 1 6Z9- 1 

with w n j_ i defined by Eq. (|17f) . As it is seen from this formula, the VEV of fermionic current is 
a periodic function of A\Li with the period of the flux quantum <I>o = 2ir/e (2irhc/e in standard 
units). It is antisymmetric about a r = 1/2. In the absence of the gauge field the VEV of 
fermionic current vanishes for special cases of untwisted and twisted fields. Of course, this result 
directly follows from the symmetry of the problem for these special cases under the reflection 
z r — > —z r . Note that L^?(0|j r |0) is a function of the ratios Li/L r and mL r . As expected, in 
the large mass limit, mL r ^> 1, the fermionic current along the direction z r is exponentially 
suppressed. 

Let us consider asymptotic limits of the VEV of fermionic current. For large values of L r , 
L r /L\ 3> 1, the main contribution comes from the term with n = 1 and to the leading order we 
have 

(0|/|0) = NLr t { ^ ] 4 P+1)/2 e~ L ^, (20) 

where w§ = Ez= P +i,^r( 27r A/^) 2 + ™ 2 and A = min(|nj + 

In the limit when the length of the one of the compactified dimensions, say z , I ^ r, is large, 
Li — > oo, the main contribution into the sum over n; in Eq. (|19j) comes from large values of |nj| 
and we can replace the summation by the integration in accordance with 



1 +°° 1 PC 

— /(27r|n, + a,|/L,) = - / 



dyf(y). (21) 



The integral over y is evaluated by using the formula from Ref. |22j and from Eq. (|19p the 
corresponding formula is obtained for the topology R p+1 x 

Now let us consider the limit when the length of one of the compact dimensions, say z D , 
is small compared with L r : Ld <C L r . In this case, in the summation over ud the main 
contribution comes from the term with minimum value of \nu + ajj\. If the parameter cud is an 
integer, the dominant contribution comes from the term with no = — old (zero mode along the 
direction z D ) and from (|19p we obtain: 

(0\f\0) = ^(0\j r \0) BPx(sl)D - l - P , (22) 

where (0\j r \0) RPx ^ s i-jD-i- P is the fermionic current in (D — l)-dimensional space with spatial 
topology R p x (S 1 ) D ~ 1 ~ P and with the lengths of compact dimensions . . . , Ad— l- In Eq. 
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I, 5d = 1 for even D and 5d = 2 for odd D. If the parameter &£> is non-integer and 
PdLt/Ld 3> 1, the argument of the modified Bessel function in Eq. (fl~9j) is large. By using the 
corresponding asymptotic formula, to the leading order we find 

/QlflON- 4iVLr " P y- sin(27ma r ) (p+1)/2 nfe 

W 1 ^ (2vr)(f+3)/2K 2^ n (p+l)/2 & nj_ 2 e ' ^ 

v ' q n=l n r q _ 2 ezi- 2 

where 

b 2 n r_ 2 = (27r(3 D L r /L D ) 2 + ^ [27r(n i + d / )L r /L i ] 2 + LV ; (24) 

Z=p+l,^r 

with rig_ 2 = (n p+ i, . . . , n r _i, n r +i, . . . , rt£)_i) . In this case the VEV of fermionic current is 
exponentially suppressed. 

In the special case with a single compact dimension we have p = D — 1, q = 1, u n r i = m, 
and the general formula f)19|) simplifies to 



<0|/|0> 



2iy m (g+i)/2 

(2tt)( D + 1 )/ 2 Z^—y—rj ( nLr )(D-l)/2 



For a massless field this expression takes the form 



jVT(CD + l)/2) ^ sin(2vrnd r ) 
<°l-? 1°) = <L 7715 • ( 26 ) 



For odd values D the series in this formula is summed in terms of the Bernoulli polynomials 
Bd(x) and one finds 

2(D-l)/2 D/2 

(VV» = (-V mi)/2 r{D/2+1)L n B»&), (27) 

for ^ a r ^ 1. In figure [J we plot the VEV of fermionic current in the simplest Kaluza-Klein 
model with D = 4 as a function of parameters a r and mL r . In Kaluza-Klein type models the 
fermionic current with the components along compact dimensions is a source of cosmological 
magnetic fields. 



3 Zeta function approach 

In this section, for the evaluation of the VEV of fermionic current we follow a different route based 
on the zeta function method [7] 123} [24"]. This allows us to obtain an alternative representation. 
To start we note that the mode sum for the fermionic current can be written as 

N + °° 

<°™ = 72- E (nr + 5r)C(l/2), (28) 

r n r =—oo 

where the generalized zeta function 
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1.0 



Figure 1: The VEV of fermionic current, L^(0|j r |0), in the model with spatial topology B? x S 1 
as a function of the phase a r and the parameter mL r . 



is introduced with the notation 



m 2 + [27r(n r + a r )/L r ]" 



(30) 



As it follows from Eq. (|28p , for the evaluation of the renormalized VEV of fermionic current we 
need to have the analytic continuation of the zeta function at the point s = 1/2. 

With this aim we first integrate over the wave vector along the uncompactified dimensions: 



D 



COO 



T(s-p/2)L r 
Uit)pI 2 T{s)V„ 



P/2- 



£ { 2 [2^(n / + d / )/L z ] 2 + 



(31) 



An exponentially convergent expression for the analytic continuation of the multiseries in Eq. 
()3ip can be obtained by using the generalized Chowla-Selberg formula [25]. The application of 
this formula to Eq. (|31|) gives the following result 



COO 



in 



D-2s-l 



r(s — {D — 1) /2) 2 1 " s m?- 2s - 1 



+ 



(4vr)( D - 1 )/ 2 T(s) (2vr)( B - 1 )/ 2 r(s) 

x 2^ cos(2vrrig_ 1 • Q g _i)/ (D _ 1)/2 _ s (m r 5(Lg„i,n^_ 1 )), 



(32) 



with L g _i = (L p+ i, . . . ,L r _i,L r+ i, . . . L D ) and = (a p+ i, . . . , a r _i, a r+ \, . . . a D ). In Eq. 

([32|) we have introduced the notation f u (x) = K u (x)/x u . The prime on the summation sign in 



321) means that the term n„ 



should be excluded from the sum and 

D 



1/2 



(33) 



i=p+l,^r 



The part in the fermionic current containing the second term on the right-hand side of 
Eq. (|32p is finite at the physical point and, hence, the analytic continuation is needed for the 
part with the first term alone. In order to do this, we apply the summation formula (I15p to 



the corresponding series over n r . After transformations similar to those already used in the 
derivation of Eq. (|19|) and by making use of the standard properties of the gamma function, 
one finds 

r(s - (D - 1) /2) 2vr(n r + a r ) /L r 



2L r (4^)( D - 1 )/ 2 r(s) — ^ ,„ 



Y2 .,,2.--+ I 



(47r) D/ 2r(g) |>(2vrna r ) ■ (34) 



As it can be easily checked, the right-hand side of this relation is finite at s = 1/2. Combining 
Eqs. (|28p . (I32p . (|34p . for the VEV of fermionic current we find the following representation 



= (27r) ( g+ l)/2 L Sm ( 27rn ^) (nLr )(D-l)/ 2 + (27r )D/2-l L ^ + ^ 
v ; n=l \ i / \ I n r =— oo 

xmf~ 2 ^ cos(27rrig_ 1 • a 9 _i)/ D/ 2_i(m r 0(L 9 _i, n^)). (35) 

Note that in the limit Li — > oo, I ^ r, the second term on the right-hand side of this formula van- 
ishes and we obtain the vacuum fermionic current in the model with a single compact dimension. 
The latter coincides with Eq. (|25j) . 

Formula (|35p is further simplified by using the relation 

+ 0O 

^2 [2vr(n r + a r )/L r ] mf ~ 2 f D /2-i( m ry) 

n r =—oo 

l~2 +°° 

= J-L 2 .m D+1 ^ n r sin(27rn r Q r )/( D+1 )/ 2 (m v / y 2 + n 2 L 2 ). (36) 

n r = l 

This relation is obtained by integrating the Poisson's resummation formula ^2n=-oo F(x)b{x — 
n ) = Tl=~oo F(x)e 2i7Tnx with the function F(x) defined by the left-hand side of Eq. fl36|). The 
integral for the right-hand side is evaluated using the formula from Ref. |22| . By taking into 
account Eq. (j36| . from Eq. (j35|) we find 



2Nm D+1 L °° 
(°b' r |°) = ( 27T )(D+i)/2 Yl n r sm(27rn r a r ) 

^ ' n r =l 

x ^ cos(27rn^_ 1 • a,-i)/(D + i)/ 2 H(L„n g )), (37) 
nj_ ie z9-i 

with the notation 

D y 2 

g(L q ,n q ) = ( £ Lfnf) . (38) 

The equivalence of two representations (119p and (I37p for the VEV of fermionic current is seen 
using the relation 

E cos(2vrn^_ 1 • a q -i)f( D+1 y 2 (rng(l, q , n q )) 

(2vr)(^ 1 )/ 2 L r ^ (D- 9 )/2+l„ , r , ,_ Q . 

= — T/^dTT— 2^ \, /(JJ-gJ/a+iCn^rWn;.!). (39) 
9 n;_ l6 z*- 1 
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The proof of this relation can be found in Appendix of Ref. [TB]. The advantage of the represen- 
tation (fT9|) . as compared with Eq. ([37|) . is that in the case of a massless field, for large values of 
ni the separate terms in the multiseries decay exponentially instead of power-law decay in Eq. 

(EZJ). 

4 Fermionic current in carbon nanotubes 

Carbon nanotubes have attracted much attention recently due to the experimental observation 
of a number of novel electronic properties. In this section we apply general results obtained 
above for the electrons in cylindrical and toroidal carbon nanotubes. A single-wall cylindrical 
nanotube is a graphene sheet rolled into a cylindrical shape. The electronic band structure of 
graphene close to the Dirac points shows a conical dispersion E(k) = vp\k\, where k is the 
momentum measured relatively to the Dirac points and vp ~ 10 8 cm/s represents the Fermi 
velocity which plays the role of a speed of light. The corresponding low-energy excitations can 
be described by a pair of two-component spinors, ipA and ips, corresponding to the two different 
triangular sublattices of the honeycomb lattice of graphene (see, for instance, [2j [3] ) . The Dirac 
equation for these spinors has the form 

(iVp^Do + iy l D t - m)i>j = 0, (40) 

where J = A,B, I = 1,2, and is defined in Eq. (|2|) with e = — |e| for electrons. To 
keep the discussion general we have included in Eq. (I40j) the mass (gap) term. The gap in 
the energy spectrum is essential in many physical application. This gap can be generated by 
a number of mechanisms (see, for example, [H [26j [271 EH])- In particular, they include the 
breaking of symmetry between two sublattices by introducing a staggered onsite energy [3] 
and the deformations of bonds in the graphene lattice [26] . Another approach is to attach a 
graphene monolayer to a substrate the interaction with which breaks the sublattice symmetry 
|27j . For metallic nanotubes we have periodic boundary conditions (a/ = 0) along the compact 
dimension and for semiconductor nanotubes, depending on the chiral vector, we have two classes 
of inequivalent boundary conditions corresponding to a/ = ±1/3. These phases have opposite 
signs for the sublattices A and B. 

The presence of the gauge field in Eq. (|40p leads to the Aharonov-Bohm effect in carbon 
nanotubes [29]. This effect manifests itself in a periodic energy gap modulation and conductance 
oscillations as a function of enclosed magnetic flux with a period of the order of the flux quantum. 
Similar oscillations arise in the VEV of fermionic current along compact dimensions. We consider 
the cases of cylindrical and toroidal nanotubes separately. 

4.1 The case D = 1 

We start with the simplest case D = 1 with a compact dimension of the length L\ = L. The 
corresponding phase in the periodicity condition we denote a\ = a. As it is seen below, this 
case can be considered as a model of a toroidal nanotube in the limit when the length of the one 
of compact dimensions is small. The corresponding effective two-dimensional Dirac-like theory 
is discussed in Refs. [30]. By summing the contributions coming from two sublattices with 
opposite signs of a, for the VEV of fermionic current one finds 

A 00 

1 > = — — y^cos(27ma)sin(27m$/$ )^i(m£n), (41) 

n=l 

where 27r$/<I>o = eA\L / (he) with <1> being the magnetic flux. The corresponding vector potential 
can be generated by the magnetic field perpendicular to the plane of torus and located inside 
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a coaxial cylinder with radius smaller than L/(2ir). For a massless case from here we have 
) = V F Sj=±-^-( ( ^/^ > o + 3&)/L, where we have defined the function 



7T ' — » 



sin(27rnx) 



71=1 



1 - 2{x}, 
2|{x}|-l, 



{x} > 
{x} < 



(42) 



The fractional part {x} on the right-hand side of this formula is defined in accordance with the 
Mathematica function FractionalPartfx]. In figure [2] we have plotted the VEV ([41 p as a function 
of the magnetic flux for different values of the parameter mL (numbers near the curves). The 
dashed lines correspond to a massless case. For the left and right panels a = and a = 1/3, 
respectively. The electric current corresponding to the VEV of the fermionic current is of order 
\e\vp/L. Note that the persistent currents in normal metal rings with this order of magnitude 
have been recently measured in Refs. [31j . 
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Figure 2: The VEV of fermionic current in D = 1 model for the periodicity conditions with 
a = (left panel) and a = 1/3 (right panel) as a function of the magnetic flux. 



4.2 Cylindrical nanotubes 

A single-wall cylindrical nanotube is a rolled-up graphene sheet in the hollow cylindrical struc- 
ture. For the case of cylindrical nanotube we have spatial topology R 1 x S 1 with the compacti- 
fied dimension of length L. The nanotube is characterized by its chiral vector Ch = (n w ,m w ), 
with n w , m w being integers determining the circumference in accordance with L = |C^| = 
a \/ n w + m w + n wi^w Here a = 2A6A is the lattice constant for graphene. A zigzag nanotube 
corresponds to the special case = (n w ,0), and a armchair nanotube corresponds to the case 
Ch = ( n w, n w)- All other cases correspond to chiral nanotubes. The electronic properties of 
carbon nanotubes can be either metallic or semiconductor-like depending on the chiral vector. 
In the case n w — m w = 3q w , q w £ Z, the nanotube will be metallic and in the case n w — m w ^ 3^ 
the nanotube will be semiconductor with an energy gap inversely proportional to the diameter. 

For the case under discussion D = 2 and the general formula for the VEV of fermionic 
current takes the form (N = 2, a p+ \ = a) 

71=1 
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with the notation a = a + e^42-L/(2-7r), and (0|j ' |0) = 0. In metallic nanotubes a = and 
for semiconductor nanotubes a = ±1/3. For graphene sheet we have two spinors that describe 
Bloch states residing on the two different sublattices. Summing the contributions from these 
sublattices and taking into account that, for these two sublattices the phases a have opposite 
signs, for the total fermionic current we find 

, )(c „) = a* £ cos(27ma) sin(2 ™^ o) I+^ (44) 

n=l 

where a = 0,1/3 for metallic and semiconductor nanotubes respectively. In these formulae 
27r ( i>/ < 3?o = eA<iLl(hc) with being the magnetic flux passing through the cross section of the 
nanotube. Note that in Eq. (|44p (and in the formulae below), we give the fermionic current 
for a given spin component. The total current is obtained multiplying by the number of spin 
components which is 2 for graphene. 

As it is seen from Eq. (|44p . in the absence of the magnetic flux the total fermionic current 
vanishes due to the cancellation of contributions from two sublattices. The magnetic flux breaks 
this symmetry and an effective current appears. However, it should be noted that, in general, 
the mass terms in the Dirac equation for separate sublattices can be different. In this case an 
effective fermionic current appears without an external magnetic field. In figure [3] we plot the 
VEV of the fermionic current for various values of the parameter mL (numbers near the curves) 
in metallic (left panel) and semiconductor (right panel) cylindrical nanotubes as a function of 
magnetic flux in units of magnetic flux quantum. 
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Figure 3: The VEV of fermionic current in metallic (left panel) and semiconductor (right panel) 
nanotubes as a function of the magnetic flux. 

4.3 Toroidal nanotubes 

A toroidal nanotube corresponds to a finite graphene sheet with the periodical boundary con- 
ditions along the transverse and longitudinal directions. This form of carbon structure was 
discovered in Refs. [5j. The carbon toroid is determined by its chiral, = (n w ,m w ), and 
translational, T = (p w ,q w ), vectors. The parameters (n w ,m w ,p w ,q w ) define the geometric 
structure and physical properties of toroidal nanotubes. For the geometry of a toroidal nan- 
otube we have the spatial topology (S 1 ) 2 with p = and q = 2 and the corresponding formulae 
for the VEV of fermionic current are directly obtained from the general results (|19p and (|37p 
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(on the persistent currents in toroidal carbon nanotubes see Refs. |32|). For a graphene sheet 
we have two sublattices with opposite signs of the phases a\. The total current is obtained by 
summing the corresponding contributions and one finds 



/ T\(tor) 2L r v F cos[27r (niai + n 2 q 2 )] 

n r =l ni=— oo v 1 1 z z/ 



^-sin[27r(ni$i + n 2 *2)/*Q], (45) 



exp (m yZ^n^+Ty 



PI) 



where r,l = 1,2, I ^ r, and 27r<3?//<l>o = eA[Li/(hc). As in the case of cylindrical nanotubes, due 
to the cancellation of contributions coming from separate sublattices, the fermionic current in 
toroidal nanotubes vanishes in the absence of the magnetic flux. An alternative representation 
is obtained by using formula (fT9l) : 



„ oo -t-oo , 

(/) ( tor) = J>n(2vrn r aW) £ ^(n, + af)f + L 2 m 2 

I 5=±n r =l n;=— oo 



xK 1 (n r {L r /L l )^[2Tr(n l + a\ S) )} 2 + L 2 m 2 ), (46) 



with r, I = 1, 2, I r, and 5j = ±a; + ^^/^o- This formula is further simplified when A\ = 0: 



(/)(*<*) = — ^ V sin(27rn r $ r /$o) cos(2 7 rn r a r .) V J[27r(n; + a,)] 2 + ^?™ 2 

' n=l n;=— oo 

xKi(n r (L r /Li)^J [27r(m + a t )] 2 + Lfm 2 ). (47) 

Note that in this case the component (0\j l \0)( tOT ^ is nonzero only for a\,a.i ^ 0. 

Let us consider the asymptotic limit of the fermionic current in toroidal nanotubes in the case 
L\ <C L 2 for a fixed value of mL 2 . First we consider the component (j 1 ) < - tor * > . For this component 
in Eq. (]46p one has 1 = 2. In the limit under consideration the main contribution to the series 
over 712 comes from large values and we can replace the summation by integration. The integral 
is evaluated explicitly and to the leading order the expression for (j l p tox ^ coincides with the 
corresponding result for cylindrical nanotubes given by (|44p (with L = L\). The behavior of the 
component (j' 2 )*- tor ^ crucially depends on whether the parameter o^p is integer or not. When this 
parameter is non- integer (for both 6 = ±) the argument of the modified Bessel function in Eq. 
(|47p (with r = 2, I = 1) is large and the dominant contribution comes from the term with n 2 = 1 
and from the term in the summation over n\ with minimum value of \n\ + o£p\. The VEV of 
fermionic current is exponentially suppressed: {j 2 p tm ' ~ exp [- x /(27r/3 1 L 2 /L 1 ) 2 + L 2 m 2 ], with 
Pi = min \n\ + af\. For nanotubes metallic along the direction z l with the length L\ and 

for A\ = one has = 0. In this case and for L\ <C L 2 the main contribution to (j 2 )^ tor ^ 
comes from the term with n\ = in Eq. (|46|) and the quantity Li(j 2 )( tor ) coincides with the 

corresponding result in D = 1 case (formula (f4"Tj) with L = L 2 and a = a 2 ). If a% 7^ and d^ 
is an integer (note that this can be satisfied for one of values 8) the main contribution comes 
from the term for which n\ + 5^ = 0. In this case only one of the sublattices contributes to 
the fermionic current. 

In figured] we plot the dependence of the fermionic current (j 2 ) ( - tor ^ in the massless case as a 
function of the ratio L1/L2 for different values of the phases {ot\, a 2 ) (numbers near the curves) 
and for ( I ) 2 / < I ) o = 0-2, 3>i = 0. In this case the component (j 1 )^ 01 ' is nonzero for ai )2 7= 



13 



only. As it was shown above and clearly seen from figure, for large values of the ratio L1/L2 the 
component (j 2 )*- tor ^ of the fermionic current tends to the corresponding quantity in a cylindrical 
nanotube with circumference L^- In the opposite limit of small values of L1/L2 the VEV tends 
to zero for semiconducting type periodicity condition along the direction z . Again, this is in 
agreement with the asymptotic analysis given before. 
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Figure 4: The VEV of the current for a massless fermionic field in a toroidal nanotube as a 
function of the ratio L\jLi for different values of the phases (01,02) and for the magnetic flux 
$ 2 /$o = 0.2, $1 = 0. 

The dependence of the fermionic current on the magnetic flux is presented in figure [5] for 
different values of the ratio L\jLi- The left and right panels correspond to toroidal nanotubes 
with phases (01,02) = (1/3,0) and (1/3,1/3), respectively. 
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Figure 5: The VEV of fermionic current in the massless case for toroidal nanotubes with 
(01,02) = (1/3,0) (left panel) and (01,02) = (1/3,1/3) (right panel) as a function of the 
magnetic flux for various values of the ratio L1/L2. 
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5 Conclusion 



We have investigated the VEV of fermionic current for a massive spinor field in the background 
of flat spacetime with spatial topology RP x S q . Along the compact dimensions the field obeys 
generic quasiperiodic boundary conditions (H|). In addition, we have assumed the presence of a 
constant gauge field. For the evaluation of the mode sum of the fermionic current two different 
approaches have been used. They give two alternative representations of the vacuum current. 
In the first approach, we apply to the mode sum the Abel-Plana type summation formula (115p , 
The renormalized VEV of fermionic current components along compact dimensions is given by 
formula (I19p . The time component and the components along the uncompactified dimensions 
vanish. The fermionic current depends on the phases in the periodicity conditions and on the 
gauge potential in the combination (|14p . It is a periodic function of the magnetic flux with 
the period of the flux quantum. In order to obtain an alternative representation of the vacuum 
current, in Section [3] we have followed the zeta function approach. An exponentially convergent 
expression for the analytic continuation of the corresponding mode-sum is obtained on the basis 
of the generalized Chowla-Selberg formula. The corresponding expression for the components 
of fermionic current along compact dimensions is given by Eq. (137p . The equivalence of two 
representations for the VEV of the fermionic current is directly seen by using the relation ()36[) . 
As a numerical example, in figured] we have depicted the dependence of the vacuum current in 
the 5-dimensional Kaluza-Klein model on the phases in the periodicity conditions and on the 
mass of the field. In this type models the fermionic current with the components along compact 
dimensions is a source of cosmological magnetic fields. 

In Section |4] we gave an application of the general results to the electrons of a graphene sheet 
rolled into cylindrical and toroidal shapes. For the description of relevant low-energy degrees of 
freedom we have followed a route based on the effective field theory treatment of graphene in 
terms of a pair of Dirac fermions. For this model we have D = 2 and the topologies R 1 x S 1 
and (S 1 ) 2 for cylindrical and toroidal nanotubes, respectively. Depending on the manner the 
cylinder is obtained from the graphene sheet, the phases in the periodicity conditions for the 
fields are equal to for metallic nanotubes and to ±1/3 for semiconductor ones. These phases 
have opposite signs for the two sublattices of the hexagonal lattice of graphene. In cylindrical 
nanotubes the total fermionic current is given by formula (|44p . In the absence of magnetic flux, 
the total fermionic current vanishes due to the cancellation of contributions from two sublattices. 
For toroidal nanotubes the two equivalent representations for the VEV of fermionic current are 
given by Eqs. (145p and (146p . As in the case of cylindrical nanotubes, due to the cancellation of 
contributions coming from separate sublattices, the fermionic current vanishes in the absence of 
magnetic flux. However, the mass terms for two sublattices can be different and in this case an 
effective fermionic current appears in the absence of the magnetic flux. 

In a way similar to that used in this paper, we can investigate the effects of non-trivial 
topology on the VEV of axial current in even dimensional spacetimes. It is well-known that in 
external electromagnetic and gravitational fields the chiral anomaly appears in the divergence 
of the axial current (for a review see [33J ) . However, in the problem under consideration these 
anomalies are absent as the spacetime is flat and the electromagnetic field tensor vanishes. 

Acknowledgments 

This work has been supported in part by the EU under the 7th Framework Program ICT- 
2007.8.1 FET Proactive 1: Nano-scale ICT devices and systems Carbon nanotube technology 
for high-speed next-generation nano-Interconnects (CATHERINE) project, Grant Agreement 
n. 216215. A.A.S. was supported by the Armenian Ministry of Education and Science Grant 



15 



No. 119. 



References 

[1] A. Linde, JCAP 10, 004 (2004). 

[2] R. Saito, G. Dresselhaus, and M. S. Dresselhaus, Physical Properties of Carbon Nanotubes 
(Imperial College Press, London, 1998); C. Dupas, P. Houdy, and M. Lahmani (Editors), 
Nanoscience: Nanotechnologies and Nanophysics (Springer, Berlin, 2007). 

[3] G.W. Semenoff, Phys. Rev. Lett. 53, 2449 (1984). 

[4] D.P. Di Vincenzo and E.J. Mele, Phys. Rev. B 29, 1685 (1984); J. Gonzalez, F. Guinea, and 
M.A.H. Vozmediano, Nucl. Phys. B 406, 771 (1993); Phys. Rev. B 63, 134421 (2001); H.-W. 
Lee and D.S. Novikov, Phys. Rev. B 68, 155402 (2003); S.G. Sharapov, V.P Gusynin, and 
H. Beck, Phys. Rev. B 69, 075104 (2004); K.S. Novoselov et al, Nature 438, 197 (2005); D. 
S. Novikov and L. S. Levitov, Phys. Rev. Lett. 96, 036402 (2006); E. Perfetto, J. Gonzalez, 
F. Guinea, S. Bellucci, and P. Onorato, Phys. Rev. B 76, 125430 (2007); A.H. Castro Neto, 
F. Guinea, N.M.R. Peres, K.S. Novoselov, and A.K. Geim, Rev. Mod. Phys. 81, 109 (2009). 

[5] H. Liu, H. Dai, J.H. Hafner, D.T. Colbert, R.E. Smalley, S.J. Tans, and C. Dekker, Nature 
385, 780 (1997); R. Martel, H.R. Shea, and P. Avouris, Nature 398, 299 (1999). 

[6] V.M. Mostepanenko and N.N. Trunov, The Casimir Effect and Its Applications (Clarendon, 
Oxford, 1997). 

[7] E. Elizalde, S.D. Odintsov, A. Romeo, A. A. Bytsenko and S. Zerbini, Zeta regularization 
techniques with applications (World Scientific, Singapore, 1994). 

[8] K.A. Milton, The Casimir Effect: Physical Manifestation of Zero-Point Energy (World 
Scientific, Singapore, 2002). 

[9] M. Bordag, G.L. Klimchitskaya, U. Mohideen, and V.M. Mostepanenko, Advances in the 
Casimir Effect (Oxford University Press, Oxford, 2009). 

[10] M.J. Duff, B.E.W. Nilsson, and C.N. Pope, Phys. Rep. 130, 1 (1986); A. A. Bytsenko, G. 
Cognola, L. Vanzo, and S. Zerbini, Phys. Rep. 266, 1 (1996). 

[11] G.L. Klimchitskaya, U. Mohidden, and V.M. Mostepanenko, Rev. Mod. Phys. 81, 1827 
(2009). 

[12] K.A. Milton, Grav. Cosmol. 9, 66 (2003); E. Elizalde, S. Nojiri, and S.D. Odintsov, Phys. 
Rev D 70, 043539 (2004); E. Elizalde, J. Phys. A 39, 6299 (2006); B. Greene and J. Levin, 
J. High Energy Phys. 0711, 096 (2007); P. Burikham, A. Chatrabhuti, P. Patcharamaneep- 
akorn, and K. Pimsamarn, J. High Energy Phys. 0807, 013 (2008). 

[13] J.S. Dowker and R. Critchley, J. Phys. A: Math. Gen. 9, 535 (1976); R. Banach and J.S. 
Dowker, J. Phys. A: Math. Gen. 12, 2545 (1979); B.S. DeWitt, C.F. Hart, and C.J. Isham, 
Physica A 96, 197 (1979); S.G. Mamayev and N.N. Trunov, Russian Phys. J. 22, 766 (1979); 
23, 551 (1980); L.H. Ford, Phys. Rev. D 21, 933 (1980); J. Ambj0rn and S. Wolfram, Ann. 
Phys. 147, 1 (1983); S.G. Mamayev and V.M. Mostepanenko, In Proceedings of the Third 
Seminar on Quantum Gravity (World Scientific, Singapore, 1985); Yu.P. Goncharov and 
A. A. Bytsenko, Phys. Lett. B 160, 385 (1985); Yu.P. Goncharov and A. A. Bytsenko, Nucl. 



16 



Phys. B 271, 726 (1986); Yu.P. Goncharov and A. A. Bytsenko, Class. Quant. Grav. 4, 555 
(1987); E. Elizalde, Z. Phys. C 44, 471 (1989); E. Ponton and E. Poppitz, JHEP 0106, 019 
(2001); H. Queiroz, J.C. da Silva, F.C. Khanna, J.M.C. Malbouisson, M. Revzen, and A.E. 
Santana, Ann. Phys. 317, 220 (2005); A. A. Saharian and M.L. Mkhitaryan, Eur. Phys. J. 
C, in press. larXiv:0911.12601 

[14] A. A. Saharian and M. R. Setare, Phys. Lett. B 659, 367 (2008); S. Bellucci and A. A. 
Saharian, Phys. Rev. D 77, 124010 (2008). 

[15] A. A. Saharian, Class. Quantum Grav. 25, 165012 (2008); E.R. Bezerra de Mello and A. A. 
Saharian, JHEP 0812, 081 (2008). 

[16] S. Bellucci and A. A. Saharian, Phys. Rev. D 79, 085019 (2009). 

[17] S. Bellucci and A. A. Saharian, Phys. Rev. D 80, 105003 (2009). 

[18] V.B. Bezerra and E.R. Bezerra de Mello, Class. Quantum Grav. 11, 457 (1994); E.R. 
Bezerra de Mello, Class. Quantum Grav. 11, 1415 (1994); L. Sriramkumar, Class. Quantum 
Grav. 18, 1015 (2001); E.R. Bezerra de Mello, larXiv:0907.4l39l Yu.A. Sitenko and N.D. 
Vlasii, Class. Quantum Grav. 26, 195009 (2009). 

[19] E. R. Bezerra de Mello and A. A. Saharian, Phys. Rev. D 78, 045021 (2008). 

[20] S.G. Mamayev, V.M. Mostepanenko, and A. A. Starobinsky, Sov. Phys. JETP 43, 823 (1976) 
[Zh. Eksp. Teor. Fiz. 70, 1577 (1976)]. 

[21] A. A. Saharian, The Generalized Abel-Plana Formula with Applications to Bessel Functions 
and Casimir Effect (Yerevan State University Publishing House, Yerevan, 2008); Preprint 
ICTP/2007/082; larXiv:0708.lT87l 

[22] A. P. Prudnikov, Yu.A. Brychkov, and O.I. Marichev, Integrals and Series (Gordon and 
Breach, New York, 1986), Vol. 2. 

[23] E. Elizalde, Ten Physical Applications of Spectral Zeta Functions, Lecture Notes in Physics 
(Springer- Verlag, Berlin, 1995); 

[24] K. Kirsten, Spectral Functions in Mathematics and Physics (CRC Press, Boca Raton, FL, 
2001). 

[25] E. Elizalde, Commun. Math. Phys. 198, 83 (1998); E. Elizalde, J. Phys. A: Math. Gen. 34, 
3025 (2001). 

[26] C. Chamon, Phys. Rev. B 62, 2806 (2000); C.-Y. Hou, C. Chamon, and C. Mudry, Phys. 
Rev. Lett. 98, 186809 (2007). 

[27] G. Giovannetti, P.A. Khomyakov, G. Brocks, P.J. Kelly, and J. van den Brink, Phys. Rev. 
B 76, 073103 (2007); SY. Zhou et al, Nature Mater. 6, 770 (2007). 

[28] G.W. Semenoff, V. Semenoff, and F. Zhou, Phys. Rev. Lett. 101, 087204 (2008). 

[29] H. Ajiki and T. Ando, J. Phys. Soc. Jpn. 62, 1255 (1993); Physica B 201, 349 (1994); A. 
Bachtold et al., Nature 397, 673 (1999); S. Zaric et al., Science 304, 1129 (2004); U.S. 
Coskun et al., Science 304, 1132 (2004); J. Cao, Q. Wang, M. Rolandi, and H. Dai, Phys. 
Rev. Lett. 93, 216803 (2004); B. Lassagne et al., Phys. Rev. Lett. 98, 176802 (2007); M.-G. 
Kang et al., Phys. Rev. B 77, 113408 (2008). 



17 



[30] K. Sasaki, Phys. Lett. A 296, 237 (2002); K. Sasaki, Phys. Rev. B 65, 155429 (2002). 

[31] H. Bluhm et al., Phys. Rev. Lett. 102, 136802 (2009); A.C. Bleszynski-Jayich et al., Science 
326, 272 (2009). 

[32] M.F. Lin and D.S. Chuu, Phys. Rev. B 57, 6731 (1998); M. Marganska and M. Szopa, Acta 
Phys. Pol. B 32, 427 (2001); S. Latil, S. Roche, and A. Rubio, Phys. Rev. B 67, 165420 
(2003); R. B. Chen et al., Carbon 42, 2837 (2004); K. Sasaki and Y. Kawazoe, Prog. Theor. 
Phys. 112, 369 (2004); Z. Zhang, J. Yuan, M. Qiu, J. Peng, and F. Xiao, J. Appl. Phys. 
99, 104311 (2006); N. Xu, J.W. Ding, H.B. Chen, and M.M. Ma, Eur. Phys. J. B 67, 71 
(2009). 

[33] S.B. Treiman, R. Jackiw, and D.J. Gross, Lectures on Current Algebra and Its Applications 
(Princeton University Press, Princeton, 1972); F. Bastianelli and P. Van Nieuwenhuizen, 
Path Integrals and Anomalies in Curved Space (Cambridge University Press, Cambridge, 
2006). 



18 



